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Abstract 

The scattering theory of electron transport allows for a compact and powerful description 
in terms oi = 1 Green functions, so-called circuit theory of quantum transport. A 
scatterer in the theory is characterized by an action, most generally a Keldysh one, 
that can be further used as a building bock of theories describing statistics of electron 
transport, superconducting correlations, time-dependent and interaction effects. The 
action is usually used in the form suitable for a two-terminal scatterer. 

Here we provide a comprehensive derivation of a more general form of the action that 
is especially suitable and convenient for general multi-terminal scatterers. The action 
is expressed as a determinant of a block of the scattering matrix obtained by projec¬ 
tion on the positive eigenvalues of the Green functions characterizing the reservoirs. We 
start with traditional Green function formalism introducing g'^ = 1 matrices and give a 
first example of multi-terminal counting statistics. Further we consider one-dimensional 
channels and discuss chiral anomaly arising in this context. Generalizing on many chan¬ 
nels and superconducting situation, we arrive at the block-determinant relation. We give 
the necessary elaborative examples reproducing basic results of counting statistics and 
super-currents in multi-terminal junctions. 
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1. Introduction 

The well-established and refined culture of theoretical description of electron trans¬ 
port in bulk solids and heterostructures was based on field-theoretical methods W and 
Keldysh Green functions Q . The pioneering works of Landauer and Buttiker [^, Q that 
unambiguously related electron transport and coherent scattering in micro-contacts have 
been regarded with suspicion: the genuine simplicity of their approach looked as a bar¬ 
baric intrusion to a sophisticated domain. It took time to appreciate the idea that the 
electron resistance is in fact scattering. Once the appreciation of this revolutionary idea 
was in place, a fast research progress has revealed many facets of the universality of 
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scattering approach and its relevance in the areas where its applicability was not at all 
obvious. The electic noise was understood in terms of scattering The sophistication 

came back when a state-of-the-art quantum calculation Q has demonstrated that the 
whole statistics of electron transport is defined by the scattering matrix. The approach 
has been applied to superconducting contacts. 

This also initiated research that combined Green function approaches with the notions 
of scattering and discrete elements giving rise to a bunch of so-called quantum circuit 
theories @ that are indispensable for accessing full counting statistics, transmission dis¬ 
tribution in complex scatterers, superconducting and spin transport in nanostructures. 
A starting point of this research was actually the paper of Buttiker and Beenakker about 
the noise in diffusive connector that seemed wrong to the author and thus motivated 
him to prove the statement on more solid grounds. Quantum circuit theories posses a 
remarkable degree of universality. A two-terminal scatterer in this approach is always 
described by an action 

‘5=iyTr 

2 ^ 


In I 1 -|- Ti- 


G1G2 + G2G1 -2 


( 1 ) 


where n labels transport channels of the scatterer, are corresponding transmission 
coefficients, while the matrices Gi ^2 characterize the states of the leads, steam from 
the Green functions and satisfy G^ 2 = 1- The matrix structure as well as the role 
of the action conforms a concrete situation from the great variety where the relation 
can be applied. In case of circuit theory of transmission distribution 0 G is a single- 
parametric 2x2 matrix and S is function to be minimized while in a theory encompassing 
interplay of Goulomb interaction and disorder 11| G may represent supersymmetric cr- 
model quantum fields, matrix structure includes time indexes and the action is a part of 
a path integral weight. 

In this article, we address the generalization of Eq. [I]to the case of multiple terminals. 
Such generalization has hardly been discussed in the literature. One of the reasons for 
this is the fact that in a (quantum) circuit theory a multi-terminal scatterer can be readily 
modelled with two-terminal ones and at least a single node connected to the leads by 
means of these two-terminal scatterers. For instance, this is a way to multi-terminal 
counting statistics. 12| However, such approach is not general. On mean-field level, it 
disregards random phase factors accumulated in the course of the scattering in the node. 
The statistics of these random phase factors can in principle be obtained if going beyond 
the mean-field level. However, this does not give an action for a concrete realization of 
these phase factors. Such action is especially important in the context of recent discovery 
of non-trivial topological phenomena in multi-terminal superconducting junctions [l^ . 

A proposal for such generalization has been made in |I4| in the context of understand¬ 
ing Fermi Edge singularity. In this paper, we provide a full and comprehensive derivation 
of this relation starting from the common textbook Green functions and extend it to the 
case of superconducting terminals. The result is expressed as a determinant of a block 
of the scattering matrix obtained by projection on the positive eigenvalues of the Green 
functions characterizing the reservoirs and is given by Eqs. I45I67I 

The structure of the article is as follows. In Section [5] we discuss the Keldysh Green 
functions, its extension to counting statistics and give a single-state multiterminal ex¬ 
ample. In Section [3] we consider a transport channel connecting two reservoirs without 
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scattering, compute Green functions and the action, recognize and heal a dangerous chi¬ 
ral anomaly. The generalization to many channels and scattering comes in the Section 
m We introduce the superconducting reservoirs in Section [S] and generalize the action 
to this case in Section m Further we elaborate on two basic examples for the block- 
determinant formula obtained. In Section [7] we derive the full counting statistics for 
multi-terminal transport of the normal electrons. In Section[S]we perform the projection 
in superconducting case, consider the ground state energy of the junction and derive a 
usefull formula for non-stationary superconducting current. We conclude in Section!^ 


2. Green functions: general 

We start our considerations with conventional definition Q of Keldysh Green func¬ 
tions in terms of averages of fermion creation-annihilation operators X being 

an element of a Hilbert space (for instance, space coordinate) 


= ('I'l'I'J) 


0 _ 1 
0 0 + 


zG{Xi,X2) 


G++ G+- 
G-+ G— 




0 + 0 
1 0 _ 


0± = 0(±(ti — ^ 2 )) 


( 2 ) 


Here Ti 2 = 'I'(Xi. 2 ) and ’’check” denotes the matrix structure in the Keldysh index 
i = ±. Let us specify to a general stationary nonequilibrium state where the density 
matrix is diagonal in the space of energy levels k and filling factor of this level is fk Since 
'k(t) = exp(—the Green function is diagonal in the levels and reads 


iGk{ti,t2) = exp(iefc(<2 - h)) 


0(^2 — H ) 1 

0 0(ti — ^ 2 ) 


- fk 


I 1 
I 1 


To get it in the energy representation, G{e) = f dte^^^G{t, 0) we note that 

^- =iGii-, luiGii =-m5{e-ek) 

J e + id - ek 

[ ^ -= -iGA-, ImG^ = in5{e - ek) 


e- iS - ek 

J = 27r6{e - ek) = i{GR - Ga) 


(3) 

(4) 

(5) 

( 6 ) 


where we have introduced advanced and retarded Green functions not depending on the 
filling factors. With those, we can represent the Green function as 


G = Gr 


-/ 1-/ 
-/ 1-/ 


+ Ga 


/-I /-I 

/ / 


(7) 


For proper description of the reservoirs we need what in old literature is called ” Green 
function in coinciding points”, G{X,X). In fact, this concept has little to do with 
geometric proximity of the points: rather, it represents a Green function ’’avegared” over 
a big number of similar states of a quasicontinuous spectrum of the same energy. Let 
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US generally define it as Ga 
such averaging reads 


= 'u^kGk, Wp being some positive weights. The result of 


Ga,v — 


1 - 2 / 

- 2 / 


2 ( 1 -/) 

1 - 2 / 


( 8 ) 


where v = — e)- For Green function in coinciding points, u is the density of 

states. Its dependence on e can be disregarded in all important cases. If the filling factors 
of the levels before the averaging depend on energy only, the filling factor /(e) in the 
above relation natually reproduces / of the levels. If not, /(e) is a weighted average of 
those and effective filling factor of this group of the states if they are used as a reservoir. 

The common Keldysh technique can be defined through the unitary evolution of the 
density matrix: 

p{t) = T ex.p{—iHt)p{—oo)T exp{iHt). 

The extended Keldysh technique Q is defined through a non-unitary evolution of the 
pseudo-density matrix with the Hamiltonians H± depending on the Keldysh index 


p{t) = T exp{—iH-t)p{—oo)T exp{iH+t) 


(9) 


We define the action S in terms of the trace of this pseudo-density matrix after its 
evolution over a big interval of time 'T, = Trp(7^). Common app lication of extended 

Keldysh technique is full counting statistics of electron transfers [l5j . In this case, H± = 
i7±y/2/, I being the operator of current to a certain reservoir and the Fourier transform 
of with respect to x gives the probabilities of transferring N electrons during the time 
interval T, 

Pjv = / 

Jq 27r 

We will be interested in variations of the action. Let us assume that the Hamiltonians 
have been changed by a little addition H± —>• Hpm + J^ab ^abiG^a^b- The corresponding 
variation of the action in the limit of small h is expressed in terms of the Green functions, 

dS = -i j J dth+(t)(4').(t)^'b(t))+ = 

= / dt(-/j:,(OG'--(i + 0,t) + h+(t)G^+(t-0,t)) 


( 10 ) 


The sign difference for ± in above equation is inconveniently annoying. Also, the 
evolution equation for G has this inconvenient sign difference, 


('-*) 


-iJ G = -7 


( 11 ) 


To avoid this nuisance, we better redefine Gnew = —GouTz departing from the conven¬ 
tional definition. With this redefinition, 


/ 1-/ 
/ 1 -/ 


G = Gr 

and the averaged function can be presented as 
Gav = -iTrvg; g = 


Ga 


1-/ /-I 
-/ / 


2/-1 2 ( 1 -/) 
2 / 1 - 2 / 


( 12 ) 


( 13 ) 










Importantly, = 1. The variation of the action becomes 


SS = -I dt + 0,t) + - 0, t)) = -Tr [hG] (14) 


where the trace in the last equality includes everything: the Hilbert space, Keldysh index 
and time. 

With this, we are ready to discuss the reservoirs. We do it in the context of a single 
level k. We take into account the transition matrix elements tkp from a state A: to a 
reservoir whose states are labelled by p. This brings self-energy to the equation for Gk 

^ = \tkp\‘'Gp = -i{r/2)g; T = 2TrJ2 \tpk\^Kep - e) (15) 

p p 


where we have averaged Gp over the states, disregarded the e dependence of T and the 
real part of E since the latter only causes an unimportant shift of Cfc. 

The Green function of the state k then reads 


Gk = 



1 

e- Cfc +i(T/2)g 


1 =1= g __ 

^ 2 e— efc±z(r/2) 


(16) 


If we take many states connected to the same reservoir, and average the Green functions 
over these states, we get Gav = —ini^g. This implies that the reservoir is reproduced and 
is completely characterized by the matrix g. 

Let us concentrate on the action. The formula for the variation can be rewritten as 


55 = -Tr [5(S)Gfc] , (17) 

from this it follows that 

5 = -Tr[lnGfc] (18) 

this makes sense: is the determinant of the matrix determining the quadratic Grass- 

man action of the creation and annihilation operators. 

Let us connect the level fc to a number of reservoirs labelled by i and consider full 
counting statistics of charge transfers between the reservoirs. Since there are several 
reservoirs, the self-energy is contributed by all of them. 


S = -^^(^,/2)5, (19) 

i 

S is a 2 X 2 traceless matrix, = —E^, the scalar E^ equals 

^ ( X! X! + 9j9]) 

\ i i>3 

5 = y'^ln((e-efc)2 + E2) (20) 

To incorporate full counting statistics, we transform g with proper counting fields xo 
g —>■ ^ so that the ^ of a given reservoirt becomes 

r 2/ - 1 2(1 - f)e^^ - 

^ [ 2/e-*x 2/ - 1 

5 



with this, 


M, + gfg^ = 2 + 4(/.(l - - 1) + /,(! - 

Theoretically, the action should vanish in the limit of y —> 0. However, this result 
may be difficult to reproduce without a tedious calculation. It is more convenient to 
compute the action allowing some freedom in the insignificant factors that bring extra 
X-independent terms and than to substract the value of the computed action at x = 0 
to correct for this. This is what we do now to arrive at 


S = T 


de 

2tt 


In 







( 21 ) 


This describes statistics of scattering of the electrons from one reservoir to another 
trought the level k with energy-dependent transmission coefficients T^-, 

TiTj 

“ (e-efe)2 + (E,r,/2)2- 

To remind, T is a, time interval at which the statistics is acquired. If T is much smaller 
than the scale at which / is changing (this is usually a temperature), this expression can 
be integrated over the energy. 


= iyr,/2y + yr,r,Mi - /,)(e*te-xd _ i) -^r,/2. 
V j j 


( 22 ) 


The /i(e) in the above expression are taken at e = ek- 

In fact, we have obtained a simplest model of a multi-terminal scattering: that gener¬ 
ated by transitions through a single resonant level. It is not the most general form of the 
scattering matrix. To consider this one, we need to turn to Green functions of transport 
channels. Those can be conventionally modelled assuming a Id spacial dependence 


3. Single one-dimensional channel 

Let us consider a Id channel with the Hamiltonian 

H = -ivd^ (23) 

sign of the velocity v is important here: at the moment, we assume it is positive. Similarly 
to what we have done for a single level, let us consider weak coupling of the channel states 
to those of a reservoir. This results in a self-energy part E = {—i/2T)g, where r is a 
typical time of escape from the channel to the reservoir, g characterizing the reservoir. 
Green function of the electrons in the channel satisfies the differential equation 

(e -I- ivdx + (i/2T)g) G{x, x') = 5{x — x') (24) 

The differential equation with respect to another coordinate reads: 

— ivdxG{x, x') + G{x, x') (e -|- {i/2T)g) = d{x — x') 
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(25) 




We can solve it in x deriving 


G(xi, X 2 ) = R(xi — x'i)G(xi, X 2 )-R ^(x 2 — X 2 ); R = exp ^-x-^—J (26) 

Since = 1, for any A 

exp (Ag) = exp(A) + exp(-A). 

Let us complicate the situation and consider a channel connected to two reservoirs. 
To this extent, let us introduce a simple x-dependence of the reservoir Green function: 
at X < 0, g = gi, and at x > 0 g = § 2 - Let us concentrate on the Green function 
in coinciding points, G(x) = G(x,x'). For a convenient normalization, it is prudent to 
rewrite it as G = —invQ = —i{l/2v)Q. Actually, Q{x,x') satisfies the equation 

(^-ie/v + dx + 0 ( 2^5 x') = 5{x - x') (27) 

We can express Q{x) in terms of Q(0) at the ’’boundary”, 

x<0: Q{x) =^e±^i^Q(0)i^-i(g(0)-5iQ(0)5i) (28) 

x<0: Q{x) =^e±^i^Q(0)i^ + i(g( 0 )- 52 Q( 0 ) 52 ) (29) 


To determine <5(0); note that for the solution not to diverge at the infinities we need 
to require 

(1 + ffi)g(0)(l - h) = 0; (1 - 52 )g( 0 )(l + § 2 ). (30) 

Besides, by reproducibility of the reservoir, we have to require Q{x) —>■ 51^2 at x —>■ =Foo. 
With this. 


g(o) = (i-5i) 


1 


ffi + ffi + h 


^ —(1 +51) = {^ + h )- 


1 


gi +92 gi+ h 


(l-ff 2 ). (31) 


Actually, g(x, x') has a discontinuity at x —>■ x'. One shows that g(x±0, x) = ±i + g(x). 
Let us compute the action. We vary the Green function of the first reservoir, 


7° 

6S = — dxTr [(5(S(x))G(x, x)] 

J —00 


1 

4tv 


dxTr [6{gi){x)Q{x) 


since Tr [< 5 ( 3 ) 5 ] = 0 for 3 ^ = 1 matrices, the integral is accumulated at distances of the 
order of 1/ur from the interface where g(0) changes into G, 

<55 = ^ /° dxTr [<5(3i)g(0)] = -^Tr [<5(3i)g(0)] (32) 
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Substituting yields 


1 


^5 = -Tr 




91 +h 


(33) 


Now we need to integrate over gi to determine the action. Eventually, this does not 
quite work. To make this explicit, let us do this for 2x2 matrices gi, 2 - This is general 
case, since gi ,2 can always be presented in the form of 2 x 2 blocks in the space of doubly 
degenerate eigenvalues of gig 2 + § 291 - We substitute g = g ■ a to obtain 


SS — SSn + SSa', 

^ c ^^9l: Us) /I I (-* w 

oSn = ^ In 1 + {gi,g 2 )); 

i9i+92r 2 

_ 1 {^ 91 , 91 , 92 ) 

“ 2 1 +(51,52) 

The variation of can be easily integrated. For general 5 , 


Sn = —Tr 

2 


In 


9i +92 


(34) 


where we have added factor of 2 to make sure = 0 if both reservoirs are identical. 

As a matter of fact, the term SSa in tbe action cannot be integrated. Indeed, if we 
set 52 II z, 5 i = (cos 0 ,sin 0 sin((), sin 0 cos(()), we obtain 

ro f.. sin 2 6 » 

SS>a — S<p n/-i , 

2(1 + cos 61) 

that is, d^deS d(j,deS. 

This signals the breakdown of the perturbation theory in E and is in fact a man¬ 
ifestation of chiral anomaly. A chiral channel between two different reservoirs is not 
well-defined, for instance, it would provide an infinite current between them. One can 
think of a less invasive configuration, for instance, 5 = 51 at a; —>■ ±00 and 5 = 52 in 
the finite but long interval of x In this case, the two boundaries between regions with 
91,2 produce the of opposite sign, so that the anomalous part cancels and the action 
is integrable. This is because the channel starts and ends up in identical reservoirs. An 
equivalent picture includes two channels of opposite velocities between two reservoirs. 
The Sa terms cancel and Sn add resulting in 


5 = Tr 


, 51+52 

In^— 


(35) 


If we substitute the reservoirs with counting field, we reproduce the known expression 
for counting statistics of the quantum point contact of ideal transparency, 

S = j ^ln(l + /i(l-/2 )e*^ + /2 (l-/i)e-*^) (36) 

This eventually sets the model in use: further we assume that the channels are 
grouped in pairs of opposite velocity, and the members of the pair share the same reser¬ 
voir. 
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4. Scattering 


Let us have many (pairs) of channels and describe scattering between those. Without 
loosing the generality we assume that all channels are at a: < 0 and the wave function 
amplitudes at x = 0 are related by scattering matrix, 

3 


We incorporate the index of the channel into the ’’check” index. To simplify the formulas, 
we assume 

Let us compute the Green function of the incoming electrons at a point x'. We define 
= Q{x' ± 0, x'). For the solution not to diverge at x —>■ —oo, it should satisfy 


{l + g)Q- =Q. 

At X > x', the Q evolves to Q'^{0,x') = R{—x')Q'^. Further, it scatters to the outgoing 
channels where the amplitudes should satisfy the condition with changed sign of gi. This 
gives 

(l-g)SR{-x')Q+ =0 

We multiply this equation by S~^ from the left to arrive at 

(1 - gs) R{-x')Q+ = 0; 33 = S-^gS. (37) 

We substitute = R{x')q^R~^ {x') to arrive at 

(1 + 9 )g~ = 0; (1 - h) g~ = O; 3 + - 3 " = 2 . (38) 


Substituting 3 ^ = 3 ± 1, we arrive at the answer similar to Eq. 1311 


9 = ( 1 - 5 )- 


1 


g + h 5 + 53 
Repeating the steps that lead to Eq. |3i 

5S = ^Tr 
2 


+ - - (1 + 5 ) = (1 + 53 ) 


1 


5 + 53 5 + 53 


we arrive at 


^(5)(1 -5)t 


1 


5 + 53 


( 1 - 53 ). (39) 


(40) 


Now we need to compute the contribution of the outgoing channels. There are sign 
differences in the equations indicating the opposite sign of the velocity. For instance, 
Eg HU becomes 


(1 + 5)5 =0; ( 1 - 34 ) 3 = 0 ; 3’^-3 = 2; 34 = 53 ^ ^ 

As a result of this, the contribution of the outgoing channels reads 


d5=iTr [%)(!+ 3 )^^ 

2 L 5 + 54 

We bring both contributions together in the following form 


*•5 = 5 '^' 


1 


{S5{g){l + g) + 5{g){l - g)S) ^ 

5-3 + ‘35 
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(41) 


(42) 


(43) 















Now it is time to integrate. The integration is less trivial than one could expect. 
The practical way is to do it in the basis where g is diagonal. The positive and negative 
eigenvalues of the matrix define a block structure: 


9 = 


1 0 
0 -1 



(44) 


The action after the integration (see Appendix for details) is expressed in terms of the 
determinant of the upper left block of S, 

S = Tr [In Sii] = In det (Sn) . (45) 

This formula can be presented in a basis-invariant form in a variety of ways, the most 
logical one is the following: 


5 = Tr 


In 



1 + g g i + g 


(46) 


where the first term is the projector on the down right block, and the second one is the 
projection of S on the upper left block. 


5. Superconducting reservoirs 

Let us generalize the Keldysh Green function of normal electrons 

= 


iGiXi,X2)=i 


-G++ G+- 
-G-+ G— 




0(^2 — ti) 

1 

0 

0(tl — ^ 2 ) 

-0(tl — ^ 2 ) 

0 

-1 

0(^2 — ti) 


(47) 


to superconducting state. We label electron operators with a = ±1, = 4', 'I'_i = 4''!' 

to introduce the Green function with Nambu indices 


zG,^(Ai,A2) = (^1.o4'2.-/3) 


SO 

1 

SO 

© 

1 

1 

0 

0(fl — ^ 2 ) 

—0(fl — ^ 2 ) 

0 

1—1 

1 

0(^2 — ti) 


(48) 


Given a superconducting Hamiltonian (where H is Hermitian and A = — A^) 

^ = E + ^'I'lAfci4/| + i4/,A*,4/t (49) 


k,l 


the time-dependence of 4/ is given by the Bogolyubov-deGennes ’’Hamiltonian” "HedG? 
.d 


^dt 


4'i 

4'_i 


= "HedG 


4'i 

4'_i 


; "HedG = 


H A 

At -H^ 


(50) 
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With this, in energy representation the Green function satisfies 


(e — "HbcIg) G =1. 


(51) 


This is not convenient since the usual potential comes to this equation in the form U(x)rjz, 
77 being the Pauli matrices in Nambu space. To circumvent this, we redefine Gnew = Gr]z. 
The redefined function satisfies 


[erjz — 'UBdo) G = i; ^BdC 


H A 
-At -H^ 


(52) 


We can equilidate the standard spin structure of A = iuyK by making unitary transfor¬ 
mation of G, 'HBdG with the matrix ((1 -|- rjz) + < 7^(1 — riz))l2 to arrive at 


"HudG 


H -iA 
iAt H 


(53) 


where H is the time-reversed H. Let us compute the Green function of a normal metal 
averaged over the states. It has a block structure in Nambu space: the upper block is 
the same as evaluated in m, while the lower block is obtained by replacing 
For a single state fc, it therefore reads (c.f. EalSTl) 


iG 


-exp(-*efc(t2 - ti)) 


— 0(^2 — ti) 
0 


In energy representation, it becomes 


-G 


down 


= Gr 


1 -/ / 
1 -/ / 


Ga 


0(ti — ^ 2 ) 


- (1 - fk) 


down/1 I \ 
k — 

-I -1 
1 I 


/ -/ 
-( 1 -/) ( 1 -/) 


; G^’^ = 


I 


e=F -I- Cfc ' 


(54) 

(55) 


(56) 


We notice that upon averaging the filling factor is taken at —e. With this, both blocks 
are presented as 


G = —iiriy 


1 + 


- 1 


F_ = - 


iTTvgN; 


2/(e)-I 2(1-/(e)) 
2/(e) I - 2/(e) 


I-2/(-e) 2/(-e) 

2(1-/(-e)) 2/(-e)-I 


(57) 

(58) 


If /(e) = 1 — /(—e) as it is the case of Fermi distribution at ji = 0,.F+ = iA = F, 
gN = VzF. 

Let us now compute the Green function of a superconductor connected to a normal 
metal reservoir. The reservoir produces the self-energy = —^l^gN-, we will assume 
l/r e, A. The Green function sought is defined by 


{er]z + iAr]^ - C - ^n) G = 1. 


(59) 


II 





















assuming real A. If /(e) = 1 — /(—e), the matrix can be diagonalized in Keldysh 
structure separately from the Nambu one, and the Green function can be presented with 
the projectors on the blocks of F, 




1 - F 


Ga', Gr,a = 


1 


(e ± i5)rj^ + Ary^, - ^ 


(60) 


In general case, the solution is more involved. Let us note that if |e| < |A|,G« = G^ and 
the Green function does not depend on /: there is no density of states for such energies 
and nothing is to be eqilibrated with the normal reservoir. In opposite case, the matrix 
E = erjz + Arjx can be presented as E = Ee, = 1. It is convenient to choose the signs 
in such a way that TT{r]ze) = vb is always positive. In this case, i^b = |e|/\/e^ — A^ is 
the normalization factor of BCS density of states. The advanced and retarded Green 
functions read 


qR,A ^ 


1 + e 


1 - 


(61) 


For averaged functions. 


qR,A _ ^ 


To compute the Green function, we project the Nambu structure in Eq. [S7]onto ± blocks 
of e: 


gN -A TrN 


lie 

—^9n 




(62) 


With this. 


6= E 

( T ,( t '=±1 


1 1 i ere 1 i (j'cFVa 

crF - ^ i id's 2 2 


(63) 


where matrices v± are obtained from V± by normalization such that = 1 For the 
standard structure not spoiled by any counting fields, v± = -E.V±/vb-, and the matrices 
v± have a structure ?? corresponding to the effective filling factors 


= /(e) 


XEXjvB 


i(l-/(-e)) 


1 T , 


/±(e) = (l-/T(-e). 


The effective filling factors account for charge imbalance in the superconductor that is 
induced provided /(e) (1 — /(—e)). 

The average G then reads: 


i-EvG = — 



(64) 


The formulas relating the Green functions and the variation of the action remain 
the same apart from 1/2 factor that comes with the trace over Nambu structure and 
compensates the artificial doubling of fermionic states in Nambu formalism. For instance, 
Eq. [TtI becomes 

J5 = _iTr[(5(E)Gfc] . (65) 

For spin-symmetric case, this factor cancels with the factor of 2 accounting for spin 
degeneracy. 
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6. Formulation for multi-terminal superconductor 


The Green functions describing the superconducting reservoir can be dealt with in 
the same way as we deal with the Keldysh functions: we just account for an extra Nambu 
index. The same holds for possible spin structure or for any other more exotic structure 
accounting for approximate degeneracy of electron states. The only essential modification 
of the approach developed in Section |4] is the Nambu index dependence of the scattering 
matrix. Indeed, the scattering matrix for normal electrons relates annihilation operators 

^out 

3 

For creation operators, this is transformed to 






= 


t,ont 


This defines the Nambu structure of the scattering matrix in conventional terms of ’’elec¬ 
tron” and ’’hole” scattering, 


c.. _ ^ Vz e 

— 2 b 


1• s® = S - 

2 ’’b’ ~ ^^3 


( 66 ) 


In the simplest case of energy and spin-independent scattering matrix, Gen¬ 

erally, s^(e) = (jy{s^{—e))'^ay. With this, we reproduce the same formulas as for the 
normal case: 

5 = ^Tr [In ^n] = ^ In det (^n) . (67) 

or, in basis-invariant form, 






( 68 ) 


As mentioned, the 1/2 factor compensates for the doubling of states in Nambu represen¬ 
tation. It cancels the spin doubling factor for spin-independent scattering. The formulas 
SSJ [03 provide an ultimately general and compact answer for statistics of electron trans¬ 
port. It is important to understand that the formulas suit perfectly the time-dependent 
scattering matrices, reservoirs and currents since time is just another index in the ’’check” 
matrix structure, and also for time-dependent counting fields. In this article, we will not 
go into details of time-dependent situation. Rather, we will elaborate two examples for 
stationary case when the matrix structure can be resolved separately at each energy. 


7. Elaborative example for normal multi-terminal statistics 


For normal case, we demonstrate the known relations for counting statistics in the 
limit of long measuring times. 

The Green function in a reservoir connected to the channel k is given by Ea ]13l As 
any matrix we can present it as an expansion in left and right eigen-vectors, ga /3 = 
'Yhi 9i'4’a'4’^■ The g has eigenvalues ±1 and without counting field the eigenvectors read 


V' 


L,+ 


1 

1 


; V’ 


i?,+ 


fk 

l-fk 




-(i-/fc) 

fk 


= 


-1 

1 


(69) 
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With the counting held Xk associated with this particular reservoir, the eigenvectors are 
changed by unitary transformation —>• exp(i ± To evaluate the action, 

we need to project onto +1 eigenvalue. Upon the projection, the element of scattering 
matrix is modihed to 

Sjk = = Sjk + (1 - . (70) 

With this, the action for full counting statistics becomes 

5 = ry ^ (in det(^)-In det (5)) , (71) 

where we have added the last term to make sure that 5 = 0 at x —>■ 0. There are many 
equivalent representations of this formula with matrices that have the same determinant. 
To give a useful one, let us transform the matrix S with unitary transform exp(ix/2) 
and multiply the result with from the left coming to 

Q = exp{-ix/2)Sexp{ix/2)S^ = 1 - / + fexp{-ix)Sexp{ix)S'^ 

We observe that det((5) = det(5)/det(5) so that 

det (Q) 

In this form, this is the same formula as given in seminal work Q : actually, Eq. 
work, that is much less popular than its two-terminal one-channel elaboration. 

For completeness, let us derive this elaboration as well. Let us consider 
scattering matrix in the basis of the channels 1,2, 


(72) 

(73) 
8 of this 
a 2 X 2 


S = 


ri ti2 
t2i r2 


its unitarity implying jnp -I- |ti 2 p = |r 2 p -I- |t 2 iP = 1, ?’iti 2 + ^' 2^12 = 0- The matrix Q 
thus reads 

n-^^\ /il^P(exp(i(x2 - Xi)) - 1) /i?'2^i2(exp(i(x2 - Xi)) - 1) 

[ /2r)'t2i(exp(-i(x2 - Xi)) - 1) hW(exp{-i{x2 - Xi)) - 1) _ 

Computing the determinant and substituting the result into [73] we arrive at the standart 
expression for the two-terminal full counting statistics that leads to binomial statistics 
of the charge transferred, 

5 = ry ^In (1 -h /i(l - / 2 )(exp(i(x 2 - Xi)) - 1) + /2(1 - /i)(exp(i(xi - X 2 )) - 1)) ■ 

(74) 


8. Elaborative example for a superconducting setup 


Let us elaborate on a similar example for a superconducting setup. For simplicity we 
assume the absence of charge imbalance in the superconducting reservoirs so that each 
reservoir is described by (c.f. Eq. isni) 


^ l + T’- 

^ = — 7 , — 9r 


l-F 


-gA\ 
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2 


(75) 








To project on the upperleft block of this matrix, we need to find its left and right eigen¬ 
vectors. Those separate in Nambu (latin index) and Keldysh structure, two independent 
vectors can be chosen as follows: 


= R^ (76) 

^2 = i>t~Ra-r9AA^ =A\ (77) 


and similar for right eigenvectors. To proceed, it is convenient to parametrize gA,R in 
the form that corresponds to zero superconducting phase yet suffuciently general one to 
account for possible complex energy dependence of these functions. The parametrization 
in terms of complex ’’phase” 9{e) is as follows: 


1 

—zsin0 1 

1 

■ -isinr 

1 

COS0 

1 i sin 9 

COS 6^* 

1 

i sin 9* 


(78) 


For pure BCS density of states, 9 is the solution of sin0 = e/|A| in upper half of the 
complex plane. At e < |A| it is real and differs by 7r/2 shift from the phase of Andreev 
scattering. @ At |e| > |A|, 9 = f snge -|- ig, where /r > 0, cosh^ = |e|/A. With this, the 
eigenvectors needed become 


=A^ 


e-^e•/2 ■ 

^ ~ 

■ ■ 

1 

gzS /2 

•\/2cos0*' 

gze/2 

y/2 cos9 


(79) 


We assign both superconducting phases (pk and counting fields Xk to the scattering 
matrix. With this, it becomes 




1 + 
2 


1 + Tz _ 

——'®+ + 



l-r]z /1 + Tz 

2 V 2 '®+ 



(80) 


where s± = exp(i(^/2±x/2)sexp(— z(^/2=fx/ 2). We project the matrix Sc onto the space 
spanned by the vectors Si ^2 going from 4x4 Nambu-Keldysh structure to a simpler 2x2 
structure. The result reads {Z = exp(*0)) 




Qii 

Qi2 


Q 21 

Q 22 

Qi2 — Qi2) Q 2 I = 021 ! Qll — Q22- 



r 2 

0 



MQM; M = 

"y /2 cos(0) 

0 

; Q = 

[S: 


-y/2cos(0*) 



fs+ + (1 - f)s. + Z-^ {fsl + (1 - f)sl) Z-^; (82) 

{-fs+{l-f) + {l-f)rs_f)p-^ + 

Z-^-frslil-f) + {l-f)slf) (83) 

(s_ - s+) - §1) ; (84) 

(1 - /) + slf + (s+(l - /) + s-/) (85) 

( 86 ) 


Since the determinant of M does not depend on y, we can skip then and concentrate on 
detQ, 

5 = i In detQ 
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(87) 






























This result provides the most complete description of the transport in a multi-terminal 
superconducting structure and is readily generalized to time- and energy dependent scat¬ 
tering matrix. Let us comprehend the structure of the answer by contemplating several 
limits. 

First of all, let us consider a limit of a non-superconducting circuit. Since Im0 —>■ -l-oo 
in this limit, we need to set Z~^ —>■ 0. This vanishes (5i2,2i and the determinant is a 
product of the determinants Qn and Q 22 - The matrix Qn in this limit approaches 

Qii fs+ + (1 — f)s-. 

This is equivalent to the matrix in Eq. 1701 The matrix Q 22 in this limit approaches 

Qii —s^(l — f) + s'^f 

which does not look the same. However, we note that Q 22 = Qn T transposition includes 
e —> —e (recall that we assume /(e) = (1 — /(—e)). Since the determinant of a matrix 
is the same as that of the transposed one, detQ 22 provides an equal contribution to the 
action. This double counting is removed by the 1/2 prefactor in the action discussed. 

Let us now derive the transport properties of the junction in the ground state. In 
this case, / = 0(—e, 1 — / = 0(e), /(I — /)=0 so that ($12 = 0 and again detQ = 
det(5iidet(322 = (detQii)^ Let us assume for simplicity that 9 is the same for all leads: 
this is the case of equal superconducting gaps in the leads. In this case, 

Qii = s_ -I- if e > 0 (88) 

Qii = s+ -I- ife > 0 (89) 


To proceed, we multiply the matrix with s_ at positive e and with s_|_ at negative e. Since 
the determinant of these matrices does not depend on counting fields, this operation does 
not change the action of interest that now becomes 

5 = ry ^ (0(e)lndet [l-be^^^sl^s!)] -b0(-e)lndet (90) 


Now we can concentrate on the eigenvalues of the unitary matrices exp(iA). The 

structure of the matrices is such that the eigenvalues form complex-conjugated pairs 
exp(z ± A). The eigenvalues exp(iA±) of the matrices are obviously the functions 

of the superconducting phases shifted by counting fields. Let us now consider the 
integral 

Although this is not immediately obvious, the integral is purely imaginary. To prove 
this, we note that 


/* (A) = 


de 

2tt 




1 -^g-2*e*g*±A 


de 

27r 




„2i9 Al 


since 0(—e) = —9*{e) and 


5k 
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= 0 


I+ r 



since the function is analytical in the upper half-plane of e and the integration 
contour can be shifted there to infinity where the integral vanishes. 

With this, we are ready to represent the action as 

S = + Xi}) - Eg{{(j>i - X,})) (91) 

Eg being a phase-dependent part of the ground state energy of the junction. This is a form 
expected for FCS in a general ground state [l^ and in particular in superconductor [I^. 
The ground state energy is expessed as a sum over the eigenvalues A({(/)i}) of the A- 
dependent part of the integral 


Eg — i I{^i) 

0<Ai<7r 


By substracting A-independent part he integral can be transformed to 


/■°° de 

— iI = E\/2= / —arg(cos(26*) — cosA) (92) 

7o 27r 

For purely BCS superconducting spectrum, the argument is tt in the energy inverval 
0 < E < Ea, where Ea is the (positive) energy of a discrete Andreev bound state 
determined from the equation 20{Ea) = A, and therefore E\ = Ea- In more complex 
situations, the discrete bound state is not formed and/or the contribution to the energy 
can come from the states of the continous spectrum. In any case, every eigenvalue A 
contibutes to the energy with a term given bv l92l 

To conclude the section, we derive a useful formula for the superconducting currents 
averaged over statistical fluctuations. This relation is obtained by singling out from the 
action the terms proportional to the first power of x^ 


I = -idiag ([s, A(1 - 2/) - (1 - 2f)Z*-‘^A'^Z*-‘^ + 2Z-‘^AQ^°^^) 


(93) 


where 

A={s + Z-A^Z-y "; = [s, f]Z*-^ + Zy, 

For a stationary case, the ’’check” structure in this formula is in the channel space, and 
it gives a contribution to the current in the channel i at an energy. It needs to be 
integrated over energy to get the stationary current. For time-dependent situation, one 
includes time in the ’’check” structure so that the formula gives the current Ii{t) while 
the energy-dependend </>, Z, f become the integral kernels in time. 


9. Conclusions 

In conclusion, we have provided a technical and comprehensive introduction to the 
Keldysh action formalism for a multi-terminal scatterer with special emphasis on super¬ 
conducting leads. We have derived a very general and compact formula and have 
elaborated on simple important examples to demonstrate the variety of its applications. 

I did this to commemorate Markus Buttiker, the pioneer of scattering approach to 
quantum transport, one of the fathers of this big, prosperous and fruitfully developing 
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research field. I admire not only his research merits: throughout 25 years of our acquain¬ 
tance I was appreciating much his daring to remain himself, to keep his own research 
style, research topics and idea sets in times where the close following of a quickly chang¬ 
ing scientific fashion seemed to be a must. He was also a charming personality and a 
good friend. 


Appendix A. Details of the variational integration 


Let us give here the details of the calculations between Eq. 031 and Eq. 
in the basis where 

9 = 

The variation of g anticommutes with g and in this basis generally reads 

9 = 


■ 1 

0 

; S = 

■ Sii 

S12 

0 

-1 


S21 

S22 


0 V 
W 0 


The inverse matrix playing important role becomes 


1 


c-l 

*^11 


0 


0 -s. 


-1 


gS + Sg 

With this, the variation of action becomes 

= Tr [^2iLE^n + '^512^^] 


We work 
(A.l) 

(A.2) 

(A.3) 

(A.4) 


Let us now determine how the blocks of S are transformed upon the variation of g. The 
matrix transforming G back to the diagonal form reads 


= 1 ± 


1 


0 V 
-W 0 


(A.5) 


Applying this to «S, S' —>■ LSL we observe that 


i5Sii — V §21 + S 12 W. 


Therefore, 


5S = Tr [dSiiSiY] 
which makes the integration straightforward. 


(A.6) 
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